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OF A SUPERSONIC WING FOR LEAST DRAG AND SUFFICIENT 
BENDING STRENGTH AT SPECIFIED FLIGHT CONDITIONS 
By James L. Amick 


SUMMARY 


This paper considers a wing-selection problem sometimes encountered 
in the preliminary design of supersonic airplanes and missiles . The • 
problem is to determine the span, section thickness ratio^, and skin thick 
ness or hollowness ratio of the wing of least drag when the plan form, 
section shape, wing lift requirement, and flight conditions are assumed 
known. The only structural requirement considered in the analysis is 
that of bending stress which is assumed to be carried entirely by the 
skin. An analytical method is presented by means of which the optimum 
wing dimensions can readily be obtained. 

An example of the application of the method to a diamond wing at 
Mach number 2.0, for a range of specified flight conditions, is presented 
From this example, for supersonic flight at low altitudes, steel wings , 
appear to have appreciably less drag than aluminum wings; whereas,- for 
high altitudes, the reverse appears to be true.. It is concluded from 
the example that wings with thickness ratios; hollowness ratios, or 
chords appreciably different from those given by the present method. may' 
have considerably higher drags. 


INTRODUCTION 


In the analysis of research data for supersonic wings and also in 
the preliminary design of supersonic airplanes and missiles, a need 
often arises -for a simple, quick method of estimating the optimum size 
and thickness of the wing for a specific application. The problem 
visualized here is one in which the thickness, hollowness, and size of 
the wing which will result in minimum drag .are to be estimated when the 
plan form and section shape have been previously selected and the flight 
conditions and the required wing lift are specified. Of the several 
structural criteria which the wing may finally be required to meet, only 
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the bending strength requirement .is considered in this .problem. The ' 
purpose of. this paper is_,tp present a simple analytical method for 
solving this preliminary problem. 

A simple criterion for selecting wing thickness has been given by 
Jones (reference l) . This criterion, which is based on a study of 
minimum spar depths of conventional aircrafifgives the maximum thick- 
ness of the wing ai^its-root as one -fifteenth the. distance from the root 
chord to ..the centroid of area of the wing panel, measured along the. - 
maximum- thickness line. Various combinations of plan form and profile 
shape may be - evaluated by comparing their maximum .lift-drag ratios, the 
thickness -chord ratio— determined by the criterion being used in each 
case. (The chord is selected so that the wing will have maximum lift- 
drag ratio at the design conditions, and the wing weight is neglected.) 

Jones' criterion is premised on the assumption that the wing is of 
the conventional thin-skinned construction) however, from considerations 
that - are given in the analysis of the present paper, it-appears that 
supersonic missiles: designed for flight at; low altitudes should. have 
nearly solid wings, which require considerably less over-all thickness 
than do thin-skinned wings, for the same strength. Jones '' criterion 
obviously could not - b e expected to apply for] .this type of wing. 

The method presented is approximate- in nature serving only as a 
rapid means for determining optimum wing dimensions, from the. drag stand - 
poirrt-for use in preliminary-design studies. In the final analysis, a 
more, rigorous approach must be employed which would consider such items 
as wing torsional strength, flutter, and divergence. 


SYMBOLS 


A 


A x 


thickness drag factor 



root cross-sectional area of, solid part of wing in streamwise 
direction - - • - 



B 


drag- rise factor 
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• 

c 

wing root chord 


w 

C D 

wing total-drag coefficient 



c D f 

wing skin-frict ion- drag coefficient 

. . - 


C L 

wing lift coefficient 



°w 

total-drag force of both wing panels 



h 

k l 

wing semi span 

ratio of section modulus of a solid wing section to c^ 

• 



ratio of moment of inertia about chord line to moment of 
inertia about neutral axis parallel to chord line, for a 
solid wing section 


A 

k 3 

ratio of distance between x-axis and x’-axis (fig. 1) to 
radius of gyration of solid wing section about chord line 



k 4 

ratio of distance between x-axis and x’-axis to distance from 
neutral axis to outermost fiber, for a solid wing section 



k 5 

ratio of distance between x-axis and neutral axis to distance 
from neutral axis to outermost fiber, for a solid wing 
section 



k 6 

t 2 

ratio of area of a solid wing section to ^c 



k 7 

ratio of volume of one wing panel to A^h 


• 

m 

hollowness ratio (ratio of maximum thickness of hollow part 
of wing cross section to over-all thickness of cross 
section) (see fig. 1) 

■ 


M 

design Mach number 


* 

m b 

bending moment about wing root chord 



n ■ 

normal load factor at design flight conditions 

— — 


n 

max 

maximum normal load factor anticipated 
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N ratio of weight of a solid wing of thickness given by equa- 

tion ( 5 ) (t] assumed to be zero) to weight of airplane minus 
wing (see equation (9) ) 

P ratio of ...thickness drag to drag due to lifting only wing weight 

for a solid wing of given chord (see equation (12)) 

q stream dynamic pressure 

r ratio of chordwise extent of hollow part of wing cross section 

to the wing chord (see fig. 1) _ 

R ratio of skin- friction drag to drag due to lifting only wing 

weight for a solid wing of minimum thickness when rj = 0 is 
assumed (see equation (20)) 

S total area of both wing panels 

t maximum thickness of wing at root 

V volume of the solid material in one wing panel 

Wf weight of airplane minus wing 

W w weight of both wing panels ; 

(1 - mr)~ : 

X = Cud 

Y ratio of moment of inertia about neutral axis of a solid 

profile, to that - of a hollow profile 

y C g distance- from root chord to spanwise center of gravity of 

wing panel. ... .. ' , 


y 


cp 


Z 


a 

6 


T) " 1 


distance -from root chord to spanwise center of pressure 

wing-root -sect ion modulus (moment - of inertia of root section 
divided by greatest distance from neutral axis) 

angle of attack 

density of wing material- .... 

ZC£ . 

y cp 


maximum allowable stress 
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ANALYSIS 

Outline of General Method 


Two of the basic objectives in the preliminary design of missile 
wings are: 

(1) The wing should produce as little drag as possible at the design 

Mach number, altitude, and load factor and with a given weight of air- 
plane minus wing. ... ' - — • 

(2) The wing should have sufficient strength that the maximum 
allowable bending stress is not exceeded in flight at the maximum 
anticipated load factor. 

At supersonic speeds, the drag of a wing increases rapidly with 
wing thickness. Also, the size of a supersonic wing is small compared 
with that of a low-speed wing of the same load carrying capacity if the 
wing is designed for flight at low or medium altitudes and at angles of 
attack near that for maximum lift-drag ratio. Achieving the minimum 
thickness ratio, therefore, is a primary consideration in designing a 
wing for least drag at supersonic speeds, while the wing weight is only 
of secondary importance (except in the case of wings designed for very 
high altitudes where low static pressures necessitate- large wing sizes). 
For these reasons, the solidity of a supersonic wing is likely to be 
much greater than that of subsonic wihgs. 

These considerations lead to certain simplifying assumptions that 
make possible a methodical selection of wing variables for the initial 
stage of preliminary design. Thus, for a supersonic wing of high 
solidity, it can be assumed’ that all the bending strength 6f the wing 
lies in the thick skin, that the wing weight, is composed entirely of 
the skin, and that no useful load is carried inside the wing. These 
assumptions permit the use of simple expressions, relating, the simple 
beam bending stress developed in the wing, the wing weight, and the wing 
drag to the thickness ratio, chord, and hollowness ratio (defined in this 
paper as the ratio of the maximum thickness of the hollow part of the 
wing section to the over-all thickness of the section, fig. 1) . 

If the assumption is made that the fuselage contributes no lift, the 
lift developed by the wing. at the design flight condition must equal the 
total airplane weight times the load factor. For a supersonic wing, 
since the wing is considered to carry none of the useful load internally, 
the total airplane weight is made up of a fixed weight (the design weight 
of the airplane minus wing) plus the variable wing weight. 

The present paper considers the problem of selecting the thickness, 
hollowness, and size of a supersonic wing so that It meets the two basic 
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objectives. The plan form, profile shape, type of skin thickness distri- 
bution, and the material of-hhe wing are. assumed .to have been previously 
selected, and the lift, drag, and spanwise center -of -pressure character- 
istics of the combination of plan. form and profile shape. are assumed to. 
be known for various thickness ratios. 

The problem is approached in." this paper by considering how the drag 
varies as the thickness ratio is varied when the hollowness ratio and " 
chord (size) are held constant. The drag of "a supersonic wing of given 
hollowness, ratio and chord can be divided into three parts: a thickness, 

drag, which increases approximately as the square of the thickness ratio, 
a drag due to lift which increases with thickness ratio (because of 
increasing wing weight), and a constant skin-friction drag. In order 
to satisfy the two basic objectives, therefore, the thickness ratio of 
the wing must be the smallest that is permitted by the second objective. 

By considering only the bending stresses, and. by using the simple beam 

formula, an expression can be obtained for the minimum thickness ratio 
as a. function of hollowness ratio and chord. - - 

Next, by making use of linear theory and the expression for minimum 
thickness ratio mentioned previously, an equation can be written which 
gives the drag coefficient of the wing as a function of. hollowness ratio 
and chord. . A differentiation of this equation with respect to hollowness 
ratio permits evaluation of the hollowness ratio that satisfies the two 
basic objectives as a function of chord. In' order to select the chord, 
the total wing drag calculated by using the previously obtained values 
of hollowness ratio and thickness ratio can be plotted against chord. By 
locating the least-drag point7 the value of chord that satisfies the two 
basic objectives and the corresponding values; of hollowness ratio and 
thickness ratio can be found. This graphical selection of the chord is 
relatively easy, whereas an analytical solution would be very difficult 
to obtain. _ ■ . " 

In case the value of chord Is dictated by other considerations, the 
thickness and hollowness of the wing which satisfy the two basic objec- 
tives can still be determined by the foregoing procedure with the last' 
step omitted. 

Development-of Method 

Selection of thickness ratio for arbitrary hollowness ratio and 
chord . - For purposes of this analysis, a wing is considered sufficiently 
strong if the maximum bending stress at the root section given by the 
simple beam formula is at all times less than' or equal to the maximum 
allowable stress. This assumption ignores additional stresses that are 
present, particularly in the case of swept wings; however, it may be 
possible to compensate to a certain' extent for these effects by adjusting 
the maximum allowable stress. (Of course, additional criteria have a 
bearing on the structural integrity of a wing; and would have to be con- 
sidered in a final design. ) 
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With, this assumption, the thickness ratio satisfying the two "basic 
objectives must be such as to give a wing-root section modulus of 


Z 


% 


max 


( 1 ) 


where Mg By is the maximum bending moment anticipated, and c a is 
the maximum allowable stress. 


The root-section modulus Z depends on the thickness ratio and 
chord of the root section and also on the hollowness of the wing. Fig- 
ure 1 shows the general type of hollow wing section that is considered 
in this paper. The profile is of such shape that multiplication of all 
ordinates by a constant factor gives a new prof ile . of the same shape but 
of different thickness ratio. (This condition is satisfied exactly by 
any profile made up of straight lines or arcs of parabolas having axes 
perpendicular to the chord; however, profiles made up of circular arcs 
approximate the condition quite well. ) The hollow space has the same 
shape as the exterior profile but is not necessarily of the same thick- 
ness ratio. The ratio of the chordwise extent of the hollow space to 
the total wing chord is denoted by r, and is considered a function of 
the hollowness ratio m. 

In order to simplify the present analysis, both m and r are 
held constant along the span. The spanwise variation of skin thickness, 
therefore, is determined by the variation of thickness ratio and chord 
along the span. For example, a rectangular wing of constant thickness 
ratio must have constant skin thickness. However, even with the restric- 
tion of constant values of m and r a rectangular wing can still be 
designed with equal bending stress at all. spanwise stations by decreasing 
the thickness ratio and the skin thickness toward the tip. 

For the general type of hollow wing section considered herein 
(fig. l), the root-section modulus can be written 


z = 


(2) 


where 


1 + 


Y = 


kg(k 3 g - l)mr * 2 (kg - l) 


2 3 k 2 

m r - k^mr + m r 


1 - ^1 - k^ + kp^mr + k^m^r 


k^, kg, k^, k^, and k^ are constants depending only on the wing- 

section shape, t is the thickness, and c, the chord, respectively, of 
the root section. For profiles symmetrical about the chord, Y becomes 


Y = 1 - m 3 r 
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The bending moment -about the root chord is caused partly by aero- - 
dynamic forces and partly by mass forces. The aerodynamic forc-es can be 
resolved into a normal force and a chordwise force, and the chordwise 
force further resolved into a drag component ; and a lift, component. .When 
the lift componentrDf the chordwise force is neglected, the bending 
moment about the root chord is " 


Mi 


■B 


y cpK + M f) n 


2 cos a 


W w 

y C g ~2 n cos a 


where y c ■ and y cg are the distances frcm the root - chord to the span- 

wise center of pressure and center of gravity of the wing panel, respec- 
tively, W w is the total weight of both wing panels, Wj*_ is the weight 

of the airplane minus wing, n is the normal load factor, and a is the 
angle of-attack. If y n _ is assumed to be independent of both the lift 

coefficient and the thickness ratio, and cos a is assumed equal to 
unity, then the maximum bending moment - is 

M %ax " I yctVme^ 1 ^ «*) <3) 

y 

where is the maximum load factor anticipated and „ q = 1 

y cp 

The value of W w in equation ( 3 ) can be expressed in terms of the 
wing geometry. For wings having the general type of cross section 
shown in figure 1, the area of the solid part of the root cross section 

is •— - • 


Ax = k 6 J c 2 (l - mr) 


where kg depends only on the profile' shape. Then the volume of the 
solid part of one wing panel ’is 

■ V = kgl^ !"! C 3(l - mr) 

where k -j is a constant depending only oh the manner of variation of 
cross-sectional area along the', span, and h - is the wing semispan. The 
total weight of both wing panels is 

^w = 2 k 5 k 7 | C 3(l - mr)5 | (b-) 

where 6 is the. density of the wing material. 
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Substituting equations (2), (3)> and (4) into equation (l) and 
solving for t/c gives the minimum thickness ratio, in terms of chord 
and hollowness ratio, as 


t. 

c 


-7T n max T l l5: 6^7 Hr " “*■) 

2k, 




+ 


'cp 


c Sc(l - mr) 

2 Va y 


'cp 


-^ W f%ax 

c 2 2k 1 a a Y 


(5) 


Selection of hollowness ratio for arbitrary chord .- The drag coef- 
ficient of a supersonic ving may be written, within the accuracy of the 
linear theory, as 


* 



( 6 ) 


where A and B are constants that depend only on the Mach number and 
the wing plan form and profile shape, and C D _^ is the skin-friction- 

drag coefficient. (Whether the linear theory gives in this same form 
the drag of a wing with spanwise variation of thickness ratio has not 
been determined. ) The lift coefficient of the wing at design flight 
conditions is 



n(w v + W f ) 
qS 


(7) 


where n and q are the load factor and dynamic pressure, respectively, 
at design flight conditions, and S is the total area of both wing 
panels. 


Wow an expression is needed for the wing weight with hollowness 
ratio and chord as the only variables. This expression can be obtained 
by combining equations (k) and ( 5 ) to give 


W. 


w 


N 2 W f T] 


(1 - mr) 2 + 
Y 


1(1 - mr) 
Y 2 


k(l - mr)^ 
N^Y 


2 


( 8 ) 
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where 


IT = 


Pep 

2i * ! 


“max 


(9) 


If t) equals zero, equation (8) "becomes 




( 10 ) 


Thus IT is the ratio of the weight of a solid wing of thickness given 
"by equation (5), calculated by assuming tj equals zero, to the weight 
of the , airplane minus wing. (For a solid wing Y equals unity. ) 

Combining equations (5) to (9) gives the drag coefficient in terms 
of hollowness ratio and chord as 


C-n — 


BT f % 2 


(1 + tj)(1 - mr) 2 N 2 ti 2 (1 - mr)^ q 


D ' afig'* 


2J C 


+ N 2 + 


PN 2 ti 2 


(l - mr)‘ 


2 + 1 1(1 - mr)^ 4(1 - mr) 2 

W 2 r] 2 Y Y y Y 2 1T 2 tj 2 Y 


tj + 


N 2 t] 2 ( 1 - mr) 2 


2Y 


' (l - mr) + 4 (1 - mr) 2 


Y 2 


TI 2 tj 2 Y 


y + c 




(11) 


where 



( 12 ) 
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The auantity P can be shcrwn to equal the ratio of the thickness ^ ra S 
to the drag due to lifting only the wing weight multiplied by the load 
factor for a solid wing of given chord. 

Differentiating equation (ll) with respect to m gives 



5m 


BWf% 2 ^r + m 


N 2 t] 2 


(13) 


where 


Qh H 


U 2 t] 2 (1 - mr)(X + l) + X 

Y 1 - mr 


ii 


(1 - mr) 2 + -ijrg + 


(1 - mr) 2 (X + 1)N 2 tj 2 + 3 x + 2 J- + (X + 2) 


(1 + i]) (l - mr) + 


- m) -.. J(1 - mr) 2 + -I 

y j y n 2. 


2 + J^ + 
"2 
n 


gfOl I gj. + 2X + (3 + ) (i . mr) 2 
Y N 2 t) 


and 


X = 


t 1 - »>g 

(r * • I> 


Equation (13) gives the slope SC^/Sm of the curve of C D against m. 
If 3cJ5m is set equal to zero, the equation defines the values of 
, sloDe is zero. For the simpler cases the value of m obtained 

hy luting ^ equal to zero can be shown to give the minimum point on 

the curve of Cp against m. 



Thus for a given chord, the hollowness ratio satisfying the two basic objectives can be found, 
since P, N, and' T) have been previously determined, and Y and r are known as functions 
of m, from 


(l - mr)^ ■ 


Jk±JL 


«V 






(l-mr) 2 +-|2 ; r+ N2T1 


2Y 


3+ T] 


H 2 !} 2 Y N 2 ^! - mr)**" (l-mr) 2 


PH = - 


] 

u 


Ht| 2 (l - mr) 2 (X+ 1) 


Y(X+ 2) H(X+2) 


iff! m-1 g l ^~i 0--»r)(3*+2) 

V jj2 2 Y(X+2) X+2 


(14) 


For T] equal to zero', this equation becomes 


PH = - 


X 


Y^(l - mr)(H(i - mr) + fr) 


(15) 


Equation (lk) is plotted in figure 2 for wing sections symmetrical about the chord-, for 
r - m and = 0.25. The type of hollow construction determined by r = m gives uniform skin 

thickness for a symmetrical double-wedge profile. 

Selection of chord. - Equations (5) and (lk) define the values of thickness and hollowness, " 
that satisfy the two basic objectives, for arbitrary chord. How in order to select the chord, 
the drag of wings having thickness and hollowness given by these equations can be computed and 
plotted as a function of chord. The minimum drag point on this plot locates the value of chord 
satisfying the two basic objectives. 
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Procedure for Use of General Method . , ■ . 

A general procedure for selecting the thickness ratio, hollowness 
ratio, and chord satisfying the two basic objectives can now be outlined 
as follows: 

(1) Select the design Mach number, altitude, maximum and design 
load factors, type of wing hollowness, weight of airplane minus wing, 

and also the wing plan form, profile shape, material, and maximum allow-. __ 
able bending stress. The lift-drag curves for wings of various thickness 
ratios and chords, having the plan form and profile, shape chosen, jmust be 
known, as well as the spanwise center-of -pressure distance. 

(2) Assume a value of wing chord. (Approximate methods presented 
subsequently can serve as a guide in the choice of chord. In case it 
may be desired to find only the thickness and hollowness ratios satis- 
fying the two basic objectives for a given value of the wing chord, that 
value may be used here. ) 

(3) Calculate W and NP from equations (9) and (12). 

(4) Read hollowness ratio m from figure 2, or a similar plot for 
the particular t] and relation between Y, r, and m considered. 

(5) Calculate W v and C L from equations ( 8 ) and ( 7 ), respectively, 

and t/c from ... 


t 

c 




(16) 


which is equivalent to equation (5). 

(6) Obtain the value of Cjj corresponding to the calculated thick- 
ness ratio and Cl, either from linear theory, higher -order theory, 

or experiment, as desired. 

(7) Calculate the wing drag from 

\ = C D<1 c 2 (IT) 

or 

(8) Repeat steps (2) to ( 7 ) for other wing chords. Plot the values 
of Dy. obtained aginst c, and determine the minimum point on the curve 
faired through these points. The values of thickness ratio, hollowness 
ratio, and chord at this point are the values satisfying the two basic 
objectives, for the given wing plan form, profile shape, and material. 
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(9) Repeat steps (l) to (8) for other plan forms, profile shapes, 
and materials, if desired. 


Approximate Methods 

In some cases it seems likely that the general method could be 
modified by the introduction of further simplifying assumptions, without 
appreciably affecting the validity of the results. In the following 
analysis three approximate methods are developed. 

Approximate method I (n = 0) .- The assumption that rj equals zero 
brings about an appreciable simplification of the general method, as has 
already been seen in equations (10) and (15). By carrying out the 
differentiation with respect to chord, instead of hollowness ratio, a 
simple equation results for the drag of a wing of arbitrary hollowness 
ratio, with thickness ratio and chord chosen to satisfy the two basic 
objectives. This equation can be differentiated with respect to hollow- 
ness ratio so that the three wing variables can be determined without 
recourse to the trial-and-error procedure of the general method. 

With the assumption that q = 0, the wing drag becomes, from equa- 
tion (11), 




Bn 2 Wf 2 


N 


PN 

Y 


N (1 - mr) 2 c 2 g(l - mr) 


N 

3 ? 






+ C 


D f q c 2 


(18) 


When this equation is differentiated with respect to c and the deriva- 
tive set equal to zero, the value of chord satisfying the two basic 
objectives, for arbitrary hollowness ratio, is obtained as 


c 2 = 



(l - mr) 2 


+ R 


(19) 


where 



2 


R = 


( 20 ) 
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(The quantity R is the ratio of the skin-friction drag to the drag 
due to lifting only the wing weight, for a solid wing of minimum thick- 
ness, when t \ = 0 is assumed.) Upon combining equations (l8) and (19) 
there results 


\ = 


Bn% f 2 




pn , 

Y 


(l - mr)‘ 


+ R + 2 


(l - mr) 

Vy 


( 21 ) 


Equation (21) gives the drag of a wing of arbitrary hollowness ratio, 
having thickness ratio and chord chosen to satisfy the two basic 
objectives. 

Differentiation of equation (21) yields 


^w 

dm 


p p dr 

Bn W ^ r + m -r— 
- DIi u dm 


PNX 


1 - mr 


♦ (1 - "’.ft * + Vr (x ♦ 2) 


(l - mr)‘ 


+ R 


When 


55s 

dm 


is set equal to zero, there are two possible roots. However, 

the root obtained by setting r + m equal to zero merely corresponds 

dm 

to a local maximum at m = 0. Therefore, the hollowness ratio satisfying 
the two basic objectives may be found, since HN and R have been previ- 
ously determined and Y and r are known as functions of m, from 


PN 


X + - 2 - (1 - mr) 
X 


1 - mr 


€ 


(l - mr)‘ 


+ R 


( 22 ) 


For the special case where r = m and the wing section is symmetrical 
about the chord line, equation (22) becomes 


PN 




( 23 ) 


Equation (23) is plotted in figure 3- 
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The procedure for using approximate method I is summarized as 
follows:- - • - - 

(1) Calculate HP and R frcm equations (9)> (12)/ and (20), with 

values of A, B, Cj^, and y^/c ctetennined either from theory or 
experiment. .. .. .. . - 

(2) Read hollowness ratio m from figure 3, or a similar plot of 
equation (22) for the particular type of wing cross section considered. 

( 3 ) Calculate the wing chord and thickness ratio from equations (19) 
and (l6)/ respectively. 

(4) Calculate the drag of the resulting wing with thickness ratio, 
hollowness ratio, and chord chosen to satisfy the two basic objectives, 
from equation (21) . 

Approximate method II (m = 0) .- In some cases the weight of a super- 
sonic wing would probably be only a small fraction of the total airplane 
weight. In such cases, the saving in drag possible through hollowing 
the wing would be negligible. With the assumption of zero hollowness 
ratio, the chord and thickness ratio that satisfy the two basic objec- 
tives become (from equations ( 19 ) and (l6), with tj = 0) 




1 

c 


^ Vmax 

^Va 


(24) 


(25) 


The drag of a wing having these values of thickness and chord is, from 
equation (21), 


Dy = 



-^(pn + 2,/m + 2) 


(2 6 ) 


Approximate method III (W v = 0) . - The analysis can be further 
simplified by assuming that the wing weight-.can be neglected entirely. 
Then equations (6), (7)/ and" (l6) combine to give the nondimens ional 
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wing drag parameter D^nW^ (which, for this approximation, is merely 

the inverse of the lift-drag ratio, since the wing weight is neglected) 
as 


jV _ ^D _ C 


4 13L a 


hW^ 


cl 


2k-, 


9- D max 

a a n 


bc l + 


(27) 


Differentiating this equation with respect to and setting the 

derivative equal to zero gives the lift coefficient at which D^/riWf is 
a minimum as ' 



(28) 


Then the drag of wings having values of thickness ratio and chord that 
satisfy the two basic objectives is found by substituting equation (28) 
into equation (27) to give 



_S_5ep a 

c 2 c < l n max 


2k 1 cr a n 


+ 



(29) 


Equation (29) describes the lifting efficiency at design flight 
conditions of the given combination of plan form and profile shape for 
all cases in which the weight of the wing can be neglected. With equa- 
tion (29) , the efficiency of wings having various combinations of plan 
form and profile shape can be compared as a function of the single 
design parameter q.nmax /o a n the Mach number, as long as the wing 

weight is negligible. Such a comparison could be of use to the designer 
in making a compromise between lifting efficiency and other requirements 
of a particular design. 

t 

It should be observed that the maximum lift-drag ratio L/D for 
wings of given bending strength occurs at a value of C-^ (see equa- 
tion (28)) smaller than that corresponding to the max imum L/D for a 
given t/c, which is ' 



* 


To see the reason for this difference, consider how the strength of a 
wing varies with chord, for a constant value of t/c. For a wing 
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supporting a given load, the bending moment^is proportional to the chord 
(when a constant plan form is assumed). The section modulus is propor- 
tional to the cube of the chord. Thus the bending Btress given by the 
simple beam formula is inversely proportional to the square of the chord, 
and therefore directly proportional to the design lift coefficient C^. 
Now as the design is made to approach the value for maximum L/D 

for a given t/c, the increase in L/D associated with an incremental 
change in approaches zero, whereas the bending stress continues ter * 

increase in proportion to C^. Thus, the last small increment in C-^ 
before the Cl for maximum L/D is reached gives an insignificant 
increase in L/D at the expense of a finite increase in the bending 
stress. Hence, in order to obtain maximum L/D for a given bending 
stress, it is. better to use a design C L smaller than that for maximum 
L/D at a given t/c and go to a smaller t/c, thus a finite— increase 
in L/D along with the finite increase in bending stress is obtained. 

As an example, the maximum L/D at Mach number 2.0 of a 10-percent -thick 
diamond wing of the type described subsequently is 4.44. However, a 
6. 5-percent-thick diamond wing operating at a smaller C-^ where the same 
bending stress id developed gives an. L/D of 5*34. The maximum L/D of 
the 6. 5-percent-thick wing is 6. 14, but the wing is not strong enough to 
go to its maximum L/D without exceeding the assumed stress. 

It is of interest to note that equation ( 29 ) can be written as 


^w_ 

nW f 


J5?(i + x ) 3 a 

n 

32k x 


^■ E1 max 

cr a n 


+ 2 



where AR is the- aspect ratio and X is the taper ratio (ratio of tip 
chord to root chord) . Now variations in AH and X are accompanied by 
relatively small changes in A, B, and y C p/h as as the leading 

and trailing edges are supersonic and if .thejtip disturbance from one 
wing does not influence the other. It appears, therefore, that a super- 
sonic wing designed for least drag should have zero taper ratio and a 
small aspect, ratio, at least for Mach numbers appreciably exceeding 
unity. 


EXAMPLE 


In order to illustrate the foregoing methods and to show the rela- 
tive importance of the various factors In a typical case, the process of 
selecting the thickness, hollowness, and size of a diamond wing to 
satisfy the two basic objectives for a range of design flight conditions 
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at Mach number 2.0 is described. The plan fozm chosen for this example 
has an aspect ratio of Ij-.O and the profile is a symmetrical double 
wedge. (See fig. 1)-. ) 

For the diamond wing considered, the geometric constants required 
for the analysis are 


k i = ^r 

kg = 1 



k 7 


1 

3 


k 

c 


1 


S_ 

c2 


1 


y cg 

c 


1 

k 


The relation r = m, which gives uniform skin thickness for the double- 
wedge profile, is used. (Of course, the skin thickness will decrease 
linearly to zero at the wing tip, since m is assumed to be constant 
along the span. ) 

The lift-drag polars of the diamond wing are given by linear theory 
as 


2 2 

c D = 2.37(1) + 0.1^5 (C L ) + 0.005 (30) 

where a skin- friction- drag coefficient of 0.005 has been arbitrarily 
assumed. (While this value cannot be correct at all the Reynolds 
numbers occurring in the example, it is believed that use of a skin- 
friction coefficient varying with Reynolds number would have little 
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effect on the comparisons presented here.) Linear theory gives the span 
wise c.enter-of -pressure distance as ■ . * 

= 0.347 

c 


In. order to illustrate the application of the general method, con- 
sider the problem of selecting the thickness, . hollowness, and size of 
the diamond wing described previously to satisfy the two basic objectives 
for the following design conditions: '* “ . . 


M =-2.0 

°a = 40,000 pounds per square inch 

= 50,000 pounds • . 

6 =0.1 pound per cubic inch (aluminum) 

q = 18.9 pounds per square inch (approx. 20,000-ft altitude) 
^ax ~ ® 

n = 1 . 


An approximate value of the chord satisfying the two basic objectives 
can be obtained from equation (24) with N set equal to zero. (which 
corresponds to neglecting the wing weight). ^ Equation (24) becomes 



and the approximate chord is found to be 15S inches, For this chord, 
the values of N and. PN are calculated from equations (9) anti (12). 
Thus 


N = 0.107 


PN = 7. 35 
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# 

From figure 2, the hollowness ratio for a chord of 15 8 inches, corre- 
sponding to. IT = 0.107 and PN = 7-35j is determined to he 0.372. 
Then substitution in equations ( 7 ), (8), and (l6) gives 


W w = 4730 pounds 


£ = 0.04l4 

c 

c L = 0.116 


From equation (3Q) the wing drag coefficient is found to be 0.0151, which 
when substituted into equation ( 17 ), gives a wing drag, of 7100 pounds. 

The foregoing procedure must be carried out for several other values 
of the chord in the neighborhood of 158 inches. From the resulting curve 
of wing drag against chord, shown to a small scale in figure 5, the opti- 
mum chord for the given conditions can be* seen to be about 158 inches. 
Thus for the given conditions, a chord of 158 inches, a hollowness ratio 
of 0.372, and a thickness-chord ratio of 0.0414 satisfy the two basic 
design objectives. 

The process of selecting thickness, hollowness, and size may be 
carried out in like manner for other design conditions. Figure 6 shows 
a plot of the drag parameter of diamond wings having values of. 

thickness ratio, hollowness ratio, and chord chosen to satisfy the two 
basic objectives calculated by the general method and also by the three 
approximate methods, for a range of design conditions. 

Figures J, 8,' and 9 present the chords, hollowness ratios, and 
thickness ratios, respectively, as functions of altitude. These fig- 
ures show, for any altitude, the chord, the hollowness ratio, and the 
thickness ratio, that satisfy the two basic objectives. 


DISCUSSION 


Figure 6 compares the drag of diamond wings at various design 
flight conditions calculated by the general method and also by the three 
approximate methods. Approximate method I, in which T) is assumed 
equal to zero, appears to give a good approximation to the general 
method, even at the higher altitudes. Approximate methods II and III 
agree reasonably well with the general method at low design altitudes, 
but diverge greatly from the general method at high altitudes. The 
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agreement of approximate methods II and in "becomes better as the air- 
plane weight is decreased or as the -ratio ofLjallawable. bending stress 
to maximum anticipated load factor is increased. 



The difference "between the values of— wing drag calculated by 
approximate methods II and III depends "on the ratio of wing weight to 
the weight of the rest of the airplane, W w/Wf* The wing weight 
increases as the design altitude of the wing is increased, since larger 
wing sizes are required at the higher altitudes. Also, decreasing 
ag/i^a x increases ..the wing weight by requiring a greater thickness ratio 

for sufficient strength. The ratio also increases with W^. 

Thus, the difference between the- two methods increases with altitude 
and is greater, .at the higher value of 
as shown in figure ,6. 


Wf and/or lower value of a & j n su 


ax> 


Figure 7 shows how the chord selected by the general and approximate 
methods increases as the design altitude is increased, so that the wing 
will: remain at an angle of attack corresponding to the efficient part of 
the lift-drag polar. However, Tn the case of the solid wing (approxi- 
mate method II ) at W f = 50,000 pounds,' the wing weight at . high altitudes 
is so great that the least drag is obtained by making the wing smaller 
than in the other cases and flying at an angle of • attack beyond the range 
for greatest lift-drag ratio. 

Values of hollowness ratio chosen to satisfy the two basic objec- 
tives are shown in figure. 8 as a function of:design altitude, "for two 
values of . W^. At low altitudes the wing size "is small, and the wing 

weight is a. small fraction of total airplane weight; therefore, a high 
solidity, which reduces the thickness drag, can be uged without greatly 
affecting the airplane weight. On the other Jhand, at high altitudes the 
wing weight is large and a large amount of- hollowness is .required. At 
the lower value of . the wing weight is a smaller fraction of. Wf, and 

consequently the hollowness ratios are smaller. 

Figure 9 shows values of the thickness ratio chosen to satisfy the 
two basic objectives. The increase iruwing chord with increasing 
altitude noted in figure 7 allows the thickness ratios of solid wings 
(approximate methods II 'and III) to be decreased at the higher altitudes 
without exceeding the limiting stress. For wings with hollowness 
selected by' the general method or approximate, method. I, the large amounts 
of hollowness required at high altitudes prevent the thickness ratio from 
being as low" as for solid wings. . 

Figure 10 compares. the drag of aluminum and steel, wings calculated 
by the general method for the same design conditions, the maximum allow- 
able stress for steel being assumed to be exactly twice that of aluminum. 
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For low and medium altitudes (where dynamic pressure and wing loading 
are high), steel wings appear to be superior to aluminum ones, especially 
if alloy steels more than twice as strong as aluminum are used; whereas, 
for high altitudes, the reverse appears to be true. (Steel wings also 
have other advantages in connection with strength at high temperatures, 
stiffness, and so forth. ) 

Figures' 5, 11, 12/ and 13’ illustrate the effect on the drag of 
diamond wings of an arbitrary choice of chord, wing loading, hollowness 
ratio, or thickness ratio, respectively. In each of. these figures one 
of the variables (chord (wing loading), hollowness ratio, or thickness 
ratio) is chosen arbitrarily, and the other two are then selected to 
satisfy the two basic objectives. Figure 5, showing arbitrary chord, 
was obtained by the general method. The other three figures were obtained 
by methods similar to the general method. These figures show that a con- 
siderable penalty in drag at design flight conditions may result if a 
value of chord, wing loading, hollowness ratio, or thickness ratio appre- 
ciably different from that given by the present method is used. 

It should be noted that the present analysis does not allow for the 
weight of internal bracing and stiffeners that would be necessary in 
cases where the hollowness ratio is large. For this reason, the curve 
in figure 6(a) calculated by the general method probably indicates 
lower drag values for high altitudes than could actually be obtained. 

In deriving the present method, the spanwise distance to the center 
of pressure y was assumed to be constant for wings of given -plan 
form and profile shape; however, in practice y C p may vary with angle 

of attack. In such a case the value of y C p corresponding to 
would have to be used. 


CONCLUSIONS 


A method suitable for use in the first stage of preliminary super- 
sonic missile wing design has been devised. The method enables. the 
analytical selection of the thickness, hollowness, and size of a super- 
sonic wing of given plan form and profile shape, having known aerodynamic 
characteristics, such that the wing will have least drag at specified 
flight conditions and yet retain sufficient bending strength. From an 
application of this method to a diamond wing at Mach number 2.0, the 
following conclusions can be drawn: 

1. The values of thickness, hollowness, and size of a supersonic, 
wing that give least drag at specified flight conditions, depend to a . 
marked extent on the particular flight conditions. . . 
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2. Over a certain range of specified jQ^ght. conditions, the general 
method of selecting thickness, hollowness, and size of supersonic wings 
for least drag presented in. this paper cah'be replaced by simpler approxi- 
mate methods with little loss in accuracy. One of - these approximate 
methods, which neglects wing weight, gives the best lift-drag ratio of a 
given combination of— plan form and profile _sh.ape.at~ a given Mach number 

as a function of a single design parameter.-;;. . . 7 

3. For supersonic flight at low altitudes, steel wings appear to 

offer a reduction in drag over that of aluminum wings; whereas, for high 
altitudes, the reverse appears to be true. -- 

4. Wings with thickness ratios, hollowness ratios, or chords 

appreciably different from those obtained by the present method may have 
considerably higher drags. 1 • . . . - 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , June 6, 1951 
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Chora, c, inches 


Figure 5.- Drag of diamond wings of arbitrary chord, having thickness 
and hollowness ratios chosen to satisfy the two basic objectives. 
General method. Wf = 50*000 pounds; q. = 18.9 pounds per square inch 

(20,000 ft altitude); — — = 5*000 pounds per square inch. Double- 

Hmax 


wedge profile; aspect ratio, 4.0; Mach number, 2.0; aluminum wings; 
r = m; C-n, = 0.005; n = 1. 
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General method ----- - 

Approximate method I (n = 0) 
Approximate method n (m = 0) 
Approximate method HI (W w = 0) 


50.000 lb 


Altitude , feet 


Dynamic pressure, q, psi 


= 5^000 pounds per square inch. 


Figure 6.- Drag of diamond wings having thickness ratios, hollowness 
ratios, and chords chosen to satisfy the two basic objectives. 
Double-wedge profile) aspect ratio, 4-.0) Mach number, 2.0) aluminum 
wings) r = m) Cp^ = 0.005) n = 1. 
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/ \ u a 

(b) = 20,000 pounds per square inch. 

^max 


Figure 6.- Concluded 
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Altitude , feet 


On 

Figure 7 .- Chord of diamond wings. — =— = 5,000 pounds per square inch; 

n max 

double-wedge profile; aspect ratio, 4.0; Mach number, 2.0; aluminum 
wings; r = m; C D = 0.005; n = 1. 
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Figure 8.- Hollowness, ratio of diamond, wings. General method and approxi- 
mate method I (tj =0). — a - = 5>000 pounds per square inch; double- 

n max 

wedge profile; aspect ratio, 4.0; Mach number, 2.0; aluminum wings; 
r = m; Cp = 0.005; n = 1. 
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a a 

Figure 9 - - Thickness ratio, of diamond wings. — — = 9,000 pounds per 

^max 

square inch; douhle-wedge profile; aspect, ratio, 4.0; Mach number, 2.0; 
aluminum wings; r = m; Cj) = 0.005; n = 1. 
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- Al uminum , 


l— Steel. aa 


10,000 psi 


Altitude, feet 


Figure 10.- Comparison, of drag of diamond wings of aluminum and of 
steel construction, having thickness ratios, hollowness ratios, and 
chords chosen to satisfy the two "basic objectives . General method. 
Wf = 50,000 pounds; double-yedge profile; aspect ratio, 4.0; Mach 
number, 2.0; r = m; C D ^. = 0.005; n = 1. 
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100,00 a 


80,000 


Altitude, ft 


60,000 


20,000 




Wing loading, 


• , lb per sq ft 


Figure 11.- Drag of diamond wings of arbitrary wing loading, having 
thickness and hollowness ratios chosen ta satisfy the two basic 

objectives, Wp = 5 0,000 pounds; — = 5,000 pounds per square inch 

n max 

doublet-wedge profile; aspect ratio, 4 * 0 ; Mach number, 2.0; aluminum 
wings; r = m; C D _ = 0.005; n = I. 
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Hollowness ratio, m 

Figure 12.- Drag of diamond wings of arbitrary hollowness ratio, having 
thickness ratios and chords chosen to satisfy the two basic objectives. 

Wf = 50,000 pounds; — - - = 5.000 pounds per "square inch; double- 

n max 

wedge profile; aspect ratio, 4.0; Mach number, 2.0; al umin um wings; 
r = m; C D = 0.005; n = 1. 
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Thickness-chord ratio, t/c, percent 


Figure 13.- Drag of diamond wings of arbitrary thickness ratio, having 
hollowness ratios and chords chosen to satisfy "the two basic objectives. 

Wf = 50; 000 pounds j. - — — ; = 5; 000 pounds per square inch; double- 

n max ' 

wedge profile; aspect ratio, 4.0; Mach number, 2.0; aluminum wings; 
r = m; C Ef = 0.005; n = 1. 
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